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We investigate the enhancement of lepton number, energy, and entropy transport resulting from 
active-sterile neutrino conversion i^e ~* Vs deep in the post-bounce supernova core followed by re- 
conversion Vs — > t^e further out, near the neutrino sphere. We explicitly take account of shock wave 
and neutrino heating modification of the active neutrino forward scattering potential which governs 
sterile neutrino production. We find that the luminosity at the neutrino sphere could be increased 
by between ~ 10 % and ~ 100 % during the crucial shock re-heating epoch if the sterile neutrino 
has a rest mass and vacuum mixing parameters in ranges which include those required for viable 
sterile neutrino dark matter. We also find sterile neutrino transport-enhanced entropy deposition 
ahead of the shock. This "pre-heating" can help melt heavy nuclei and thereby reduce the nuclear 
photo-dissociation burden on the shock. Both neutrino luminosity enhancement and pre-heating 
could increase the likelihood of a successful core collapse supernova explosion. 

PACS numbers: 14.60.Pq,95.35.+d,97.60.Bw,98.80.-k 



I. INTRODUCTION 

There are many purely astrophysical and nuclear 
physics uncertainties in the core collapse supernova prob- 
lem. However, the weak interaction in general and neu- 
trino physics in particular play pivotal roles in nearly 
every aspect of the collapse of the core of a massive star 
and likely in any subsequent supernova explosion as well. 
It is sobering to contemplate that collapsing stellar cores 
will pass through regimes of matter density and neutrino 
flux which have never been probed in the laboratory and 
which could be affected significantly by new physics in 
the weakly interacting sector. Moreover, the existence 
of neutrino rest masses, unexplained and unpredicted by 
the Standard Model of particle physics, points directly 
at the possibility of new neutrino physics. 

In this paper we explore the effects of plausible exten- 
sions of the Standard Model in the weakly interacting 
sector on models for the explosion mechanism for core 
collapse supernovae. In particular, we investigate the ef- 
fects of an electroweak singlet ("sterile") neutrino Vs on 
the physics of energy and lepton number transport in the 
supernova core and on the process of shock re-heating. 
The ranges of sterile neutrino rest mass and active-sterile 
vacuum mixing angle investigated here include those 
parameters of interest for sterile neutrino dark matter 
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and pulsar kicks and related issues [l^, [20|, l2l|. The 
LSND experiment [22,,^,^] and recent mini-BooNE ex- 
periment [2^ do not constrain the sterile neutrino mass 
and mixing parameters considered in this paper. 

The general features of core collapse supernova evolu- 
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tion are dictated largely by entropy considerations [26l |. 
Stars with initial masses in excess of ^ 10 Mq evolve 
quickly to their evolutionary endpoint: a low entropy core 
supported by relativistically-degenerate electrons and, 
therefore, subject to dynamical instability. The collapse 
of this core is halted at or just beyond the point where nu- 
clear density is reached. The gravitational binding energy 
released in this prompt collapse and in subsequent quasi- 
static contraction is more or less efficiently converted into 
seas of neutrinos of all kinds. The "bounce" of the core 
generates a shock wave which moves out. However, the 
energy in this shock is sapped by the photo-dissociation 
of nuclei passing through it. This process is an inevitable 
consequence of the substantial entropy jump across the 
shock front and of basic nuclear physics. 

The details of the mechanism or mechanisms whereby 
the deleterious effects of nuclear photo-dissociation are 
ameliorated, a viable shock is re-born, and an explosion 
originates remain elusive. However, ever since the work 
of Bethe and Wilson [27! the broad outlines of a solution 
are plausibly clear. The prodigious energy in the neu- 
trino and antineutrino reservoirs in the collapsed core is 
radiated from the surface of the proto-neutron star (the 
neutrino sphere) and is deposited in material behind the 
stalled bounce-shock, "re-heating" it and thereby driving 
a Type II, lb, or Ic supernova explosion. 

However, one-dimensional simulations of this process, 
though containing detailed treatments of the nuclear 
equation of state and neutrino transport, nevertheless 
are challenged in producing convincing explosions. Much 
recent attention has focussed on multi-dimensional hy- 
drodynamic, convective, or acoustic enhancement of neu- 
trino energy transport [2^ [2^, [s^, HH above the neutrino 
sphere as a means of augmenting neutrino heating of mat- 
ter below the shock. These schemes succeed in producing 
explosions. However, as yet they do not include the level 
of sophistication in, e.g., neutrino transport and nuclear 
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equation of state employed in the one-dimensional mod- 
els for all relevant regimes of time and space. 

Our previous work [s^ on the effects of active-sterile- 
active {ve — *■ Vs — * ^e) neutrino flavor transformation 
in the in-fall epoch of supernova core collapse suggested 
a means by which neutrino energy transport could be 
augmented. Conceivably, this could be a solution to the 
shock re-heating problem. However, a key uncertainty 
not addressed in Ref. [32| was the effect on this process 
of the shock wave itself. Here we will tackle this issue. 

In Section II we summarize the salient features 
of active-sterile-active neutrino flavor transformation 
physics and its effects during the in-fall epoch. In Sec- 
tion III we consider the ways in which the shock wave 
modifies the thermodynamic conditions which help deter- 
mine how sterile neutrino production and re-conversion 
proceed. We also discuss sterile neutrino induced "pre- 
heating" and the possibility of a reduced nuclear photo- 
dissociation burden on the shock. In Section IV we dis- 
cuss shock re-heating and the enhanced prospects for 
a supernova explosion which could be a by-product of 
active-sterile-active neutrino conversion schemes. We 
give conclusions in Section V. 



II. IN-FALL PHASE NEUTRINO FLAVOR 
CONVERSION 

In this section we briefly summarize our previous work 
[3^ on the effects of active-sterile neutrino flavor conver- 
sion on the in-fall phase of a core collapse supernova. 
The key result of this earlier work was the discovery 
that electron neutrino conversion into a sterile neutrino 
species Ue — s- Vs could feed back on electron capture 
(e~ + p ^ n + Ve) during collapse and alter the potential 
governing flavor transformation so as to produce a double 
Mikeyev-Smirnov-Wolfenstein (MSW) resonance 
It is this double resonance structure which can lead to 
the the re-conversion of the sterile neutrinos. With such 
a double resonance arrangement, at least some electron 
neutrinos will experience Vf, ^ Ug ^ as they move 
from higher toward lower density in the core. 

For simplicity, we consider 2x2 neutrino flavor mixing 
where, in vacuum, we have 



= cos6'|i^i) + sin6'|i^2), 
Ws) = - sin 6* I i^i) -I- cos 6*11^2) 



(1) 
(2) 



Here 9 is an effective 2x2 vacuum mixing angle for the 
i>e ^ i^s channel, and \i>i) and \i'2) are light and heavy, re- 
spectively, neutrino energy (mass) eigenstates with mass 
eigenvalues mi and rn2, respectively. The relevant mass- 



squared difference is 6m 



Since we will be 



concerned with sterile neutrino rest mass scales ~ keV, 
we will have m2 S> mi, and so Sm^ « m| = m^. 

An electron neutrino (i^e) propagating coherently in 
the medium of the core will experience a potential 
stemming from forward scattering on all particles (elec- 
trons/positrons, nucleons/quarks, and other neutrinos) 



that carry weak charge. This potential is 

o /T\ / 1 /I 



(3) 



where rib = pNj^ is the baryon number density, p is the 
density in gcm~^ and Na is Avogadro's number, Gf is 
the Fermi constant, and the net lepton abundances rela- 
tive to baryons are, e.g., Ye = (Ue- — ne+)/nh with, e.g., 
rie- the electron number density. The terms proportional 
to Y^^ , Yi, , and Y^^ in this potential stem from neutrino- 
neutrino forward scattering and must be corrected for the 
non-isotropic nature of the neutrino distribution func- 
tions at locations which are above the neutrino sphere 
[35| . At any point inside the star or above it, electron 
antineutrinos, i.e., f'e's, will experience a potential with 
the same magnitude as that experienced by i^e's, but with 
opposite sign. Of course, sterile neutrinos experience 
no forward scattering potential. 




FIG. 1: Right panel shows the core density profile with radius 
r, while the corresponding profiles for MSW resonance energy 
-Eros (solid) and Ve chemical potential fi^^ (dashed) are shown 
in the left panel. Here i?res takes its minimum value -B^s" at 
vq. For a particular neutrino energy, an MSW resonance can 
occur at two locations (densities), e.g., ri (pi) and r2 (P2). 



At a given location, a neutrino (i^e or Vs) with en- 
ergy Ej-cs will experience an MSW medium-enhanced res- 
onance where 



6m^ cos 29 
2V 



2V 



(4) 



Physically, this is the neutrino energy where the effec- 
tive in-medium mass associated with the active neutrino 
matches the rest mass associated with the sterile state, 
ms. The last approximation in Eq. (4) follows for the 
reasons given above and because the vacuum mixing 
angles we consider here are very small {e.g., satisfying 
sin^ 29 ~ 10"^). 

In medium the forward scattering potential will mod- 
ify not only the effective masses of the active neutrinos 
but also the unitary relation between the neutrino flavor 
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states (weak interaction eigenstates) and the (instanta- 
neous) mass eigenstates |j^i(i)) and |t'2(i)), where t rep- 
resents any Affine parameter along the neutrino's world 
line. We can express the in-medium transformation in 
direct analogy to that in vacuum, 

\iye) = cos9M{t)\iyiit)) +smeMit)W2{t)), (5) 
= -sineM{t)Wi{t))+coseM{t)W2{t)). (6) 

A similar unitary transformation applies to the antineu- 
trinos but with a different mixing angle 0m (0- 
active-sterile neutrino oscillation scenario where neutrino 
transformation is enhanced and antineutrino transforma- 
tion is suppressed, at resonance we will have 0m(^ics) = 
7r/4, i.e., maximal mixing. The region in space where the 
effective in-medium mixing angles 9u (or 0m) are large 
and near maximal is termed the resonance width. This 
width is 5t « jdlnF/dip^ tan 26* and so is expected to 
be small for the neutrino parameters and conditions we 
treat here. 

So long as the the neutrino mean free paths are large 
compared to the MSW resonance width, we can regard 
neutrino flavor evolution as coherent, at least as far as 
the application of the MSW formalism is concerned [3^ . 
This is true even when the active neutrinos are trapped 
and thermalized in the core. Note, however, that at very 
high densities, such as those we expect to encounter deep 
in the core near and after core bounce, this condition will 
break down. There may be so many scattering targets 
for the active neutrinos in this case that the neutrino 
mean free paths are comparable to or shorter than the 
MSW resonance widths. We term this the incoherent 
or scattering-dominated case. In this regime, scattering- 
induced de-coherence of the neutrino fields will dominate 
the conversion of neutrino flavors. In particular, this can 
be the case for the Ve i^s channel of most interest 
here. Note, however, that since the de-coherent neutrino 
(antineutrino) flavor conversion rate is proportional to 
sin^ 26*14 (i) (sin^ 2^M(i)), the potential V and the MSW 
resonance condition still play a significant role in deter- 
mining the locations where this conversion is significant. 
Ref. [4)] and references therein discuss this physics in de- 
tail, while Ref. [s^, [s^ discusses uncertainties and con- 
troversies associated with de-coherence in high density 
matter. 

Employing a simple nuclear liquid drop model [2^ . [38l | 
and degenerate electron equation of state in a one-zone 
homologous collapse code [33|, we found the double res- 
onance structure discussed above. Fig. [1] gives a graphic 
summary of these results. The equation of state and one- 
zone collapse code employed in obtaining these results is 
discussed in the Appendix of Ref. ^32] . 

These calculations also showed that near the surface of 
the core, where the density is p ^ lO^^g/cm^, the MSW 
resonance energy Et-cs for v^, ^ Vs tends to be much larger 
than the v^. chemical potential (Fermi energy) ^.^^ . Pro- 
gressing inward from the edge of the collapsing core, Ey-^s 
first decreases while increases continuously. Near the 



density p ~ 10^'^g/cm'^, E^cs and /i^^ become compara- 
ble and large-scale v^. — > Vg conversion starts. Once this 
conversion process begins in earnest, i?rcs increases with 
further increases in density. In this latter phase of the col- 
lapse, Et-cs stays slightly above , with both quantities 
increasing with increasing density. As will be discussed 
in the next section, a feedback process keeps E-^ca hover- 
ing just above . Ultimately, at core bounce, when the 
collapse is halted, the matter density is near nuclear mat- 
ter density {p ~ IQ^'^ g/cm'^) and the relevant neutrino 
energies are large since p^^ ~ 150 MeV. Fig. [1] illustrates 
these trends. 

Our earlier work jsl] speculated that the double MSW 
resonance structure could facilitate enhanced neutrino 
energy, entropy, and lepton number transport from deep 
in the core to regions nearer the proto-neutron star sur- 
face {i.e., the neutrino sphere). Essentially this enhance- 
ment comes about because a neutrino, initially a v,, in 
our case, will spend part of its time as a sterile neutrino. 
While it is in the sterile state, this neutrino will move 
at almost the speed of light. As a result, the effective 
mean free paths and diffusion coefhcients for these neu- 
trinos will be re-normalized upward. Interestingly, our 
estimates suggested that the best prospects for transport 
enhancement through this mechanism could be obtained 
with sterile neutrino mass and vacuum flavor mixing pa- 
rameters which overlap the ranges of these that give vi- 
able sterile neutrino dark matter P, 0, S, 0, S H 13, B B 
M ini, [H, E Q [B d, 113, d • 

However, a significant caveat on these conclusions is 
that the calculations of Ref. [s^l dealt only with the in- 
fall epoch of core evolution. The bounce shock generated 
near the edge of the homologous core could be expected 
to move outward, through the outer core, and modify the 
thermodynamic variables and composition in this region. 
These modifications, in turn, could be expected to alter 
the forward scattering potential which governs sterile 
neutrino production and/or re-conversion. 



III. EFFECT OF SHOCK WAVE PASSAGE 

Assessing the impact of post-shock active-sterile-active 
neutrino flavor transformation requires adroit attention 
to a few key issues in supernova shock formation and 
propagation. As the initial iron core collapses, an inner, 
homologous core will maintain a roughly self-similar, in- 
dex 3 polytropic structure (26l . [39j . This makes intuitive 
sense because the pressure support in the star is domi- 
nated by relativistically degenerate electrons with Fermi 

1/3 

level (chemical potential) /ie ~ 11.1 MeV(pio Ye) , 
where pio is the density in units of 10^*^ gcm^^. 

However, as electron capture proceeds and the pres- 
sure is relatively reduced, only a smaller, "inner core" 
can continue to collapse in this self similar and homolo- 
gous manner. Homology (in-fall velocity proportional to 
radius) allows a one-zone calculation to be meaningful, 
as each location in the inner core will experience a por- 
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tion of a common temperature, density, and composition 
history ,26] . 

The remainder of the initial iron core which is above 
and outside the inner core is termed the "outer core." The 
inner core is essentially an instantaneous Chandrasekhar 
mass Afic ~ (^e)^- When the central density reaches 
the point where nucleons touch (nuclear density), this 
core will bounce as a unit and serve as a piston. The 
shock will form at the edge of this inner core. The ini- 
tial shock energy will be of order the gravitational bind- 
ing energy of the inner core and will scale as (Yg)^'^^'^ 
[38|. As a result, there is some uncertainty in this ini- 
tial shock strength depending on nuclear and sub-nuclear 
density equation of state, composition, and electron cap- 
ture physics issues. In broad brush, however, we expect 
the entropy-per-baryon S (in units of Boltzmann's con- 
stant /cb) to jump by a few units at the shock front. 

This entropy jump can be significant because the core's 
material during the collapse itself, as well as the un- 
shocked material in the outer core ahead of the shock, is 
characterized by low entropy, S ~ I. In the lower density 
regions of the outer core, an entropy jump AS* > 3, for 
example, is usually enough to shift the nuclear composi- 
tion in Nuclear Statistical Equilibrium (NSE) from heavy 
nuclei to free nucleons and alpha particles. We will re- 
fer to this phenomenon as nuclear photo-dissociation or 
nuclear "melting." 

As the shock propagates through the outer core and 
melts nuclei it loses energy. This is because each nu- 
cleon is bound in a nucleus by '--^ 8 MeV. This represents 
10^^ ergs (= 1 Bethe) per 0.1 Mq of material transiting 
the shock front. Since the shock is born with an energy 
^ 1 Bethe and the outer core mass may be ~ 0.7 Mq, nu- 
clear photo-dissociation quickly degrades the shock into 
a "dead," standing accretion shock. 

Whether subsequently the shock can be re-energized 
by, e.g., direct or convectively- or hydrodynamically- 
enhanced neutrino heating or electromagnetic or acoustic 
energy transport remains an ope n q uestion as discussed 
in the Introduction [H, [H, [30, [3l[. However, by any 
objective standard, the energy (~ 1 Bethe) in observed 
Type II supernova shocks/explosions is small compared 
to the energy (~ 10 Bethe) in the neutrino seas initially 
trapped in the core, and miniscule compared to the en- 
ergy 100 Bethe) in the neutrino seas a few seconds 
post-core-bounce. Active-sterile neutrino transformation 
can tap into this reservoir and change the way in which 
neutrino energy is transported in and around the super- 
nova core. 

As discussed in the last section, direct active-sterile- 
active neutrino flavor transformation could re-normalize 
upward the neutrino energy transport rate, thereby in- 
creasing the neutrino luminosity at the neutrino sphere 
and so boosting the shock re-heating rate. Also, the ef- 
ficacy of the various re-heating schemes may depend on 
how far out the shock progresses before it stalls. In turn, 
this depends, among other variables, on electron cap- 
ture and the shock energy remaining after nuclear photo- 



dissociation in the outer core. (See the discussion on this 
point in Ref. [13] ■) Any effect like pre-heating which 
diminishes the nuclear photo-dissociation burden could 
translate into a larger stall radius for the shock, in turn, 
helping to increase the effectiveness of the various shock 
re-heating processes. 



A. Feedback between resonance energy and Vs 
Fermi level 

An important finding in the calculations of Ref. [s^l 
was that the active-sterile MSW resonance energy E^es 
exhibited a minimum which was located well inside the 
core. The density profile and i?rcs profile at bounce is il- 
lustrated in Fig. [TJ The location of the minimum in Ey-cs 
at bounce is another way to divide the core. As a con- 
sequence of this minimum, the first resonance Ve t^s 
may occur in the inner core, while the re-conversion res- 
onance, the second one, I's — > i^e, typically occurs in the 
outer core. Note that at the inner resonance, inside of 
the location of the minimum in -Ercs, the i^e Fermi en- 

1/3 

ergy fi^,^ w 11.1 MeV(2pio Y„J ' tracks just below i^rcs, 
increasing with increasing density just as does i?rcs- 

Another key finding of Ref. [32] was that in the region 
inside of the resonance energy minimum there is a feed- 
back between sterile neutrino production, Yi,^ , and Y^. 
which keeps -Bros tracking just above /i^^. This feedback 
process is a result of the high degeneracy in the electron 
neutrino distribution function. If the system were per- 
turbed so that Ej-es were lower than fii,^ , there would be 
prodigious sterile neutrino production which would tend 
to lower the local net electron lepton number and return 
the system to a state with Bj-cs > Mi/^ • 



B. Shock wave modification of sterile neutrino 
production 

The passage of the shock through a region can alter 
the relation between Ej-cs and and so can infiuence 
sterile neutrino production there. As long as i^ros stays 
well above the electron neutrino Fermi energy ^i,^ , the 
production of sterile neutrinos is negligible. However, 
the shock wave can supply heat/entropy and can cause a 
discontinuous change of physical quantities {e.g., density 
and entropy). Immediately behind the shock front, we 
might expect the density jump to result in a smaller gap 
between Ej^^s and /i^^ . This could be accompanied by 
enhanced production. However, as outlined above, we 
expect this condition to be temporary, as the feedback 
effect will push Bics above again. 

To take into acount this effect in our one-zone calcula- 
tion, we added heat and entropy "by hand." Specifically, 
to simulate the conditions in newly shocked regions of 
the core, we instantaneously increased the density by 
Ap = 10^^ g/cm^ and the entropy-per-baryon (in units 
of /cb) in three different cases by AS* ^ 0.6, 2, 3 as 
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measured at density p = 10^^ g/cm'^. We assume that f3- 
equilibrium and Nuclear Statistical Equilibrium (NSE) 
are attained instantaneously. This will be a decent ap- 
proximation in the very high density regions where the 
first MSW resonance will be located, e.g., inside or just 
outside the inner core. 

The entropy increments AS" that we employ are cho- 
sen to be values characteristic of the early stages of 
Shockwave formation. These values are smaller than the 
AS* ~ 10 entropy jump across the shock which is ex- 
pected at later times or larger radius. However, our val- 
ues make sense in a rough, physical sense: For a Chan- 
drasekhar mass initial iron core (~ 10^^ baryons) col- 
lapsing to nuclear saturation density, we expect an in-fall 
kinetic energy at bounce ~ 10^^ erg which, if dissipated 
as heat at temperature T ~ 1 MeV, would give AS" 1 . 
(See the discussions in Ref. [26| and Ref. |38j.) Going 
beyond this crude estimate is tricky. 

As best we can ascertain, our values of AS* at relevant 
locations and epochs in the core bracket the results of 
some published large-scale and detailed numerical simu- 
lations. Both Ref. JlJ and Ref. [12] seem to infer val- 
ues of AS* for relevant locations and epochs which are 
within the range we consider here. However, as we will 
see below, within this range of entropy jump there can 
be significant differences in i/g — > i^^ ^ i^e effects. 

We calculate i/g production and the infiuence of this 
process on the core in the following manner. First, we 
prepare an initial density profile. This is meant to be 
characteristic of the core just prior to core bounce. We 
take this profile to be that of a self-similarly contracted 
(homologous) index n = 3 polytrope with central density 
Pccntrai = 3 X 10^'*g/cm'^. We then choose the location 
of the shock front on this profile and take the density 
there as the initial density when the shock front arrives. 
We take the other initial physical quantities from the 
results of our in-fall one-zone calculation at the initial 
density. We then apply our increments in density and 
entropy. Following a numerical procedure similar to that 
used to get the initial model, we use the results of an 
appropriate one-zone calculation to get the new, post- 
shock thermodynamic and lepton number quantities for 
the given increments Ap and AS". We use these altered 
conditions to estimate the production of sterile neutrinos 
and the feedback of this process on the potential V. 



C. Heating of the outer core 

For a given neutrino energy, we can identify the lo- 
cation of the second, outer resonance by using one-zone 
collapse calculation results for the run of potential V (or, 
equivalently, Sros) and the corresponding density pro- 
file. In order to assess the effects of neutrino flavor re- 
conversion — > Ve in this outer region, we need to esti- 
mate how many i/g's are delivered and how much energy 
is deposited at the second resonance. 

This can be estimated by assuming adiabatic neutrino 



flavor evolution through MSW resonances. (Ref. [S^l dis- 
cusses why adiabatic evolution is a good approximation 
here.) In the adiabatic limit we can assume that all I'e's 
contained in neutrino energy range AEj^cs, corresponding 
to the MSW resonance potential width SV, are converted 
to sterile neutrinos 1/^- The width of the resonance in ra- 
dial coordinate is Sr = (dr/dV) SV — Titan 20. Here the 
potential ("density") scale height is Ti. — \d\nV/dr\~^ . 
Another expression for the spatial resonance width is 
Sr = HAE^cs/ Ej-es- Making use of the resonance con- 
dition Eq. we can express this as 



Sr 



2V^AE„ 



dV 



dr 



(7) 



Using this, we can show that the re-conversion rate per 
baryon for i^g ^ i^e at the second resonance is related to 
the corresponding rate per baryon for Ve i^s conversion 
at the first resonance by 



^istPist dV/ dr \2nd 
P2nd dV/dr list 



(8) 



' 2nd 



At any location we can designate L = Yg + Y,^^ as the to- 
tal electron lepton number per baryon. Neutrino flavor 
conversion — ^ Vg (vg v^.) produces a negative (pos- 
itive) time rate of change of this quantity, L — dL/dt, 
respectively. In employing Eq. ([5]), we evaluate dV/dr 
numerically using the in-fall one-zone calculation profile. 
In this equation, pist (p2nd) and rist (7'2nd) are the density 
and the location of the first (second) resonance, respec- 
tively, as illustrated schematically in Fig. [TJ The energy 
transfer rate per baryon from the first to the second res- 
onance obeys a relationship in obvious analogy to that 
in Eq. dSD. 

At the location of the second, outer resonance we take 
account of the heat and lepton number deposited by 
Vs Vf, by re-running the one-zone code with these up- 
dated quantities but with the density fixed at its original 
value. This gives us estimates of the change in thermody- 
namic variables that accompany this "pre-heating." We 
continue this calculation of z/g — > I's — > f^e energy transfer 
to locations in the outer core until the shock wave reaches 
the position tq where E^^s takes its minimum value (see 
Fig. [ID. 

Fig. O shows the profiles of entropy, temperature, and 
heavy nucleus mass fraction ATh at the completion of 
this energy transfer process. Three profiles are shown, 
corresponding to three different shock strength scenarios 
with entropy jump (as measured at p = 3 x 10"'^^ g cm~'^) 
AS" = 0.6 (triangles), AS" = 2 (squares), and AS" = 3 
(circles), respectively. We may view these profiles as 
snapshots of conditions when the shock front is located 
at tq. The figure also includes the results of the origi- 
nal in-fall calculation for comparison. Fig. [5] shows that, 
depending on shock strength and the £^res profile, ster- 
ile neutrino-induced pre-heating could result in at least 
partial 50%) melting of heavy nuclei in the outer re- 
gions of the core ahead of the shock. This could represent 
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FIG. 2: Effects of the shock in the interior (curves on the 
right) and shock-modified Ve. Us ^ Ue. pre-heating in the 
outer core (curves on the left). Heavy nucleus mass fraction 
Xh, temperature k^T (in MeV), and entropy per baryon S (in 
units of fee) are shown for three different cases. These cases 
correspond to three different shock strength scenarios with 
entropy jump (as measured at p = 3 x 10^^ g cm~^) AS — 0.6 
(triangles), AS = 2 (squares), and AS = 3 (circles), respec- 
tively. Results of the pre-shock, in-fall one-zone calculation 
are also included (dotted lines). The minimum of the reso- 
nance energy is located around p = 7 x lO^^gcm""^. This 
location divides the curves on the left and right, as described 
in the text. 



IV. SHOCK RE-HEATING 

In their pioneering work on core collapse supernovae, 
Mayle and Wilson [4^ Q obtained vigorous explosions 
in the late-time shock re-heating model, even in one di- 
mension. This result was, and continues to be, at odds 
with the results of other detailed one-dimensional sim- 
ulations, some more sophisticated in their treatments 
of the nuclear equation of state and neutrino transport 
[2i,[33,liSliEl43,H|4i,[5S[p[5l[53. Mayle and 
Wilson got their result by invoking neutrino convective 
transport in the core to increase the neutrino luminos- 
ity at the neutrino sphere. Though the physical basis 
for this effect {i.e., their "neutron fingers") has been re- 
pudiated, their result taught us a valuable lesson: The 
efficacy of neutrino heating in rc-cnegizing the stalled 
shock is a sensitive function of neutrino and antincutrino 
transport in the core and the corresponding luminosities 
at the neutrino sphere. The process of z/g ^ v^. 

flavor conversion in the core could be just the sort of 
neutrino energy transport augmentation that could aid 
the core collapse supernova explosion process i32l] . 



A. De-coherent production of sterile neutrinos 
inside the proto-neutron star 

Neutrino flavor evolution deep in the central region of 
the post-bounce core will be coUisionally-dominted. The 
characteristic density in the central core at this epoch 
will be near or above nuclear saturation density, p ^ 
3 X 10^^ g cm"-^, and scattering- induced de-coherence will 
be the primary channel through which sterile neutrinos 
are produced from the seas of active neutrinos [13] . 

The total (left-handed plus right-handed Dg) ster- 
ile neutrino emissivity £ (energy emission per unit mass 
per unit time) can be estimated by employing average 
neutrino and antincutrino flavor conversion probabilities 
{Praise ^ '^s]P,t)) and {Pm{ve ^ Us;p,t)), respectively, 
as functions of neutrino or antincutrino momentum p and 
location parameter t, energy-dependent neutrino and an- 
tincutrino scattering cross sections (in principle on all 
weakly interacting targets) au^{E) and ap^(E), respec- 
tively, and integratin g o ver neutrino and antineutrino 
fluxes and energies E |M . [ssl . [H, , 



a substantial reduction in the nuclear photo-dissociation 
burden for the shock. Even though our estimates are 
schematic in nature and crude on a quantitative level, 
this result is sufficiently dramatic that it is clear that 
the existence of sterile neutrinos in the mass and mix- 
ing ranges discussed here could alter the the energetics 
of core collapse supernova shock propagation. 



£ K — ( d<^>,^Ea,^{E)hp„,{iye ^ i^s;p,t)) 
rriN J 2 

+ — f d^p^Ea^SE)\{Pm{ve^i>s]P,t)), (9) 
rriN J 2 

where mAr is an atomic mass unit (essentially, the average 
free nucleon mass). In the conditions of near weak and 
near thermal equilibrium in the post-bounce central core, 
the differential neutrino and antineutrino fluxes d<^^^ and 
c?$p_. (or number densities duu^ and drip^), respectively. 
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can be expressed as 



1 



(27r)3e^/^-=-"-= +1 
1 E'^dE 



(10) 



d<i>. 



cdn,-. 



d^p 



1 



(27r)3ei5/T,,-^„, 

1 £;2^^; 



(11) 



where the (f'e) degeneracy parameter is r/^,^ = ^i^^/T^^ 
(jype — [i-VclT^v^ ~ ~^i^e)i respectively. The neutrino and 
antineutrino temperatures T^^ and Tp^, respectively, are 
essentially the same as the matter temperature. Here the 
speed of light is c. The average oscihation (transforma- 
tion) probabihties in Eq. ^ are given by 

1 A{Ef sin^ 26 



2 A{E)^ sin^ 26* + £|2 + [A{E) cos 26* - F]^ ' 



(12) 



{Praise Vs]P, t)) 

1 A(£:)2 sin^ 26* 



2 A{EY sin^ 261 + _^ [A(£;) cos 26 + V]'^ 

(13) 

FoUowing Ref. [H^l, and for the purpose of simple esti- 
mation, here we will take the Ve and scattering cross 
sections to be those appropriate for free nucleons. These 
are roughly 

o-^. {E) w CTp^ {E) w IMGlE^. (14) 
In Eqs. and we employ the notation 

A(]5) EE (5mV2p« m2/2£;« A(£;). (15) 
The quantum damping rate for neutrinos is 



D = V,J2 = J d^,^a,^{E)/2. 



(16) 



The analogous quantum damping rate for antineutrinos, 
D, has a form directly analogous to that for D. 

The effect of the de-coherent Vs and i?s production on 
the potential V has been studied in the context of a col- 
lapsed stellar core in Ref. [s^l- There it was argued that 
V should evolve toward zero on a time scale short com- 
pared to the characteristic proto-neutron star core dy- 
namical time scale. Accordingly, we shall take = in 
Eq. (fT2)) and Eq. (fTS]) in the following discussion. This 
will facilitate a simple estimate of the sterile neutrino 
emissivity deep in the central region of the proto-neutron 
star after bounce. 



Enhancement of neutrino luminosity behind the 
shock 



We have estimated the effects of shock passage on ther- 
modynamic and composition variables in the outer parts 
of the core by employing one-zone simulations of shock 
propagation through these regions. In doing this, we use 
the same numerical procedure described in Section III for 
gauging the effects of shock passage in the inner parts of 
the core. However, in the case of the outer core, we take 
account of the v^. i/g pre-heating of the material 

prior to the arrival of the shock. Therefore, our initial 
conditions for shock passage in the outer core for this 
calculation are chosen to be those given by the Vs i^e 
energy deposition process described in Section HI and 
shown in Fig. [51 

The results are intriguing. For the case of a strong 
initial shock (AS" = 3 as measured at density p = 
lO^'^ gcm~3)^ our calculations show that the double reso- 
nance structure characteristic of the in-fall regime is de- 
stroyed. In this case, however, the resonance energy E-^cs 
remains well above the i^e Fermi energy fi^,^ . This, in 
turn, suggests that any i^e which is converted to a sterile 
neutrino Vs by scattering-induced de-coherence deep in- 
side the core, yet possesses an energy above the value of 
Ercs at the neutrino sphere, will encounter an MSW reso- 
nance further out, nearer the neutrino sphere, and will be 
coherently and adiabatically re-converted to a ly^ there. 
Fig. [3] shows the results of the one-zone calculations that 
suggest this scenario. 

It can be seen in Fig. [3] that both the E-^^s and fi,^^ 
curves are monotonic with increasing density and each 
has positive slope. Therefore, the highest energy neutri- 
nos will tend to deposit their energy {i.e., be re-converted 
to J'e's) deepest in the core. This could result in more 
heating by i^s ^ i^e transport enhancement with in- 
creasing depth which could, in turn, promote convective 
instability and further augmentation of neutrino energy 
transport. 

In any case, since our estimates show that the res- 
onance energy E^es asymptotes out to about E^^^e ~ 
100 MeV at the outer edge of the core, we can conclude 
that the Ve's converted to sterile species in the inner re- 
gions of the core where /^^^ > E°^^c reconverted 
to i^e's prior to escaping the core. On account of the 
quadratic energy dependence of the Vg absorption cross 
sections, such re-converted high energy i^eS are certain to 
deposit their energy and be thermalized on times scales 
short compared to any transport time scale. 

Our calculations suggest that a weaker initial shock 
will not eliminate the double resonance structure left at 
the end of the in-fall epoch. Fig. 3] is analogous to Fig. [3] 
but shows the results of a one-zone calculation with ini- 
tial shock strength AS* = 0.6. In this case neither the 
i^e i^s pre-heating of the outer core or the shock 

passage event itself can change composition, density, and 
temperature enough to disrupt the general form of the 
runs for E^^^ and [i^^ . 
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p [g/cm3] 

FIG. 3: One-zone calculation results for resonance energy 
Eres (in MeV) and Ve chemical potential (Fermi energy, 
in MeV) are shown as functions of density p (in gcm""^). Cir- 
cles and squares represent E^cs and ^v^, respectively. Filled 
symbols correspond to the values of these quantities for an 
assumed strong shock {/S.S = 3, as described in Section III). 
This case includes the effect of — > i^e reconversion and asso- 
ciated pre-heating ahead of the shock, as well as the effect of 
the shock itself. The effect of post-bounce pre-heating alone 
is shown by the quantities with the open circles and squares. 
For comparison, iJres (dashed line) and (dotted line) are 
given for the in-fall (pre-shock, no pre-heating) case. 




1013 10' 

p [g/cm3] 

FIG. 4: Same as Fig. 3, but for the case of a weak shock 
(AS' = 0.6). 



We conclude that there may be a threshold in shock 
strength beyond which the double resonance structure 
at the end of in-fall is replaced by the single outer reso- 
nance regime in Fig. [3] What is this threshold in shock 
strength? The answer to this question is hard to get at 
with our simplistic model. However, a fair guess based 
on our one-zone scheme with its liquid drop equation of 
state would be AS > 2 (as measured at p = 10^'^ g cm"'^). 

This is significant but, ultimately, unsatisfying because 
large-scale numerical supernova simulations, depending 
on the initial model and on in-fall physics, may pro- 
duce initial bounce shocks with strengths below, near, or 
above this threshold. For example, the calculations by 
the Mezzacappa group [i^l , appear to produce shocks 
with strengths AS* « 2 by our measure. This would be 
near or above the threshold for erasing the in-fall epoch 
double resonance structure. However, the simulations by 
the Janka group [131 suggests a range of shock strengths 
which could be near the threshold. This issue has to 
be resolved before we can be confident of the effects of 
Ve ^ i^s ^ Ve on core collapse supernovae. 

Note, however, that in either the weak or strong shock 
case, sterile neutrinos produced at high energies deep in 
the core could be converted to z^e's further out. This is all 
we need to enhance energy deposition behind the shock 
and, therefore, increase the shock re-heating rate. All 
that remains is an estimate of this heating rate. This re- 
quires an estimate of the sterile neutrino cmissivity deep 
in the core. 

Following Ref. [s^l, we can get a rough estimate of 
the energy radiated in sterile neutrinos Vg per unit mass 
and per unit time - the emissivity - in the region of the 
core where the v^, Fermi energies are /j,^^ > 100 MeV. As 
outlined above, we take y = deep inside the proto- 
neutron star and approximate the Fermi distribution as 
a step function, i.e., completely degenerate, with degen- 
eracy parameter 77^^ ^ 1. In this limit, flavor conversion 
in the channel v^. — > Vg gives rise to sterile neutrino Vs 
emissivity 



1.66 



■ sin^ 29 



dE 



1 + iD^E'^/mt 



1.66 G?, . 

sin^ 29 , 

IGn^ruN V 2D 



(17) 



where we ignore contributions to the emissivity stem- 
ming from Pg's. Noting that the integration parameter 
^ satisfies ^ < 1 for 100 eV < < 1 MeV and that 
typically /i^^ ~ 150 MeV, we can calculate the emissivity 
to leading order in ^ to find 



f « (2 X lO^^^erg s-^g-^) 



sin^ 20. 



(18) 



This is then the rate per gram at which energy in sterile 
neutrinos is flowing out of the inner parts of the core. 
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On account of adiabatic MSW resonant Vg flavor 
conversion, the fraction of the deep core's Vg energy flux 
which is carried by neutrinos with energies above the 
resonance energy at the outer edge of the core, E'^^^'^ , 
will be deposited in the regions just below the neutrino 
sphere. Using a calculation in obvious analogy to that in 
Eq. P?|) , we can estimate the effective emissivity for this 
"re-captured" sterile neutrino energy, 

£[Vs Ve) 

xi(P™K ^ i^e;p,t)) (19) 

where, as argued above, E^^^'^ « 100 McV. Using the 
same approximations made in evaluating Eq. (jl7p . we 
find 

£i,,s - >^e) « 1.4 X 10^2 g-1 (^^^) . (20) 

Since the inner part of the core which generates the ster- 
ile neutrinos has a mass ~ IMq, the energy deposited 
per unit time near the edge of the neutron star could 
be prodigious. Of course, this conclusion depends on a 
host of active-sterile neutrino mass/mixing matrix issues 
including, e.g., the effective 2x2 angle 9 characterizing 
Ue ^ I's vacuum mixing. 

If we take nis ^ 1 kcV and sin^ 29 ~ 10^^, correspond- 
ing to the "sweet spot" for sterile neutrino dark matter 
and beneficial supernova effects picked out in Ref. [33 |. 
then the emissivity in Eq. (|20)) suggests that we could 
possibly double the j/g energy resident just below the neu- 
trino sphere. Though this energy would be deposited in 
the form of i^e's, rapid re-establishment of beta equilib- 
rium would imply that this energy is shared among all 
six active neutrino species. This energy sharing roughly 
will be weighted by the relative numbers of active neu- 
trino species in equilibrium. However, if the extra I'e's 
are deposited quite close to neutrino sphere, energy re- 
distribution becomes a difficult neutrino transport issue. 
Since there is a preponderance of fe's, we can guess that 
there will not be equal amounts of energy in the Ve, i^e, 
v^, D^, Vt, and Vt seas. 

On the other hand, since shock re-heating is mostly 
effected through the charged current capture processes 

-'r n ^ p + e~ and -\- p ^ n -\- e'^ , it is the v^. and 
luminosities at the neutrino sphere which are most 
important. We could be conservative and assume equal 
energy sharing so that the v,, and z/g seas get a third of 
the extra energy deposited by Vg — *■ near the neutrino 
sphere. In this case, and for a range of mixing angles 
relevant for sterile neutrino dark matter, we could ex- 
pect roughly a '-^ 10% to ~ 100% increase in the sum of 
the Vf, and luminosities. This, in turn, could lead to 
comparable increases in the re-heating rate of the shock. 



V. LEPTON NUMBER TRANSPORT AND THE 
ROLE OF Ai- AND r- FLAVOR NEUTRINOS 

In this section we discuss the active-sterile-active neu- 
trino flavor transformation-induced flows of electron, 
muon, and tau lepton numbers and the effects of these on 
supernova physics. The v^, — s- i/g process outlined 
above will transport electron lepton number from deep 
in the core to the vicinity of the neutrino sphere. In the 
course of describing this process, we made no consider- 
ation for mu (z/^, v^) and tau {vt, T>t) flavor neutrinos. 
Surely, if electron neutrino fiavors mix in vacuum with a 
sterile species, likely so will mu and tau fiavor neutrinos. 

In broad brush, the lepton number transport rate for 
Vf, ^ Vs ^ Vf, should dominate over the rate for Vf, — > 
Vs — > Df. and, for that matter, the rates for t/^ — s- i^s — > I'p, 
Vf_, ^ Vs ^ v^, Vr ^ Vs ^ Vr, and Vr ^ Vs ^ v-r as 
well. The argument to support this assertion is based 
on the relative populations of the various active neutrino 
species. 

Keep in mind that the inner core, the "piston" for 
shock generation at bounce, though experiencing an in- 
crease in entropy stemming from the dissipation of in- 
fall kinetic energy, nevertheless remains relatively low in 
entropy and full of its original electron lepton number 
excess. Immediately after bounce, the temperature in 
the core is T ~ lOMeV, while the i^e Fermi energy is 
Mi'e ~ 100 MeV. (See the discussion in Ref. ^5^.) In 
the standard stellar collapse model, these conditions will 
persist for of order a neutrino diffusion time scale, i.e., 
seconds. This is a time comparable to or longer than the 
shock re-heating time of interest here. 

The i^e's will have a negative chemical potential 
{—[iu^). The mu and tau flavor neutrinos must be pair 
produced, and as a consequence they will have zero chem- 
ical potential. In the conditions of beta equilibrium in the 
inner core, the number density of z^e's will be '--^ /z^ , while 
the number density of i/g's will be ~ T'^ exp (— /i^^/T), 
and the number densities of all mu and tau flavor neu- 
trino species will be ^ T^. Clearly, there should be 
a large excess of z/g's over the other neutrino species 
in the time frame of interest. As a result, during this 
time, de-coherence associated with the scattering of ac- 
tive neutrino species will produce far more sterile neu- 
trinos (i^s's) than the opposite handedness "anti" -sterile 
neutrinos (i^s's). 

The picture we have of the supernova core in this time 
frame is then as follows. We have an inner core "source" 
producing a large flux of very high energy Vs 's and lower 
fluxes of lower energy i^s's. The z^^'s will be preferen- 
tially transformed to I'e's via Vs — > z^e near the neutrino 
sphere. This is because in this region the forward scatter- 
ing potential for mu or tau neutrino conversion to sterile 
neutrinos will be negative. With a negative potential, 
only antineutrinos can be matter-enhanced. 

For example, the forward scattering potential for the 
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flavor conversion channel 



+4K 



is given by 

- 1 + 2r,, 

). 



(21) 



An analogous expression holds for the potential, V^^ , rel- 
evant for Vs ^ i^T, but with the coefflcients of Y^^^ and 
in Eq. (|2T|l swapped. Since we expect « 0.35, 
Yu^ ~ 0.05, and — Y^^ ~ initially, we will have 
Vv < and V^^ < 0. This, in turn, implies that 
only the channels Ds ^ f^/^ and Vs ^ S't, respectively, 
can be matter-enhanced and resonant near the neutrino 
sphere. These processes will be sub-dominant compared 
to Vs ^ Ve because the energy and fluxes for the 's will 
be lower than those for j/g's as argued above. 

The dominant i^g — > v^. conversion process will lead to 
the region near the neutrino sphere being "charged up" 
with positive electron lepton number. Given the energy 
emissivities discussed in the last section, for example we 
might expect an additional electron lepton number per 
baryon AY^^ ~ 10^^ ergs^V (lOOMeV • lO" baryons) w 
0.1 to be deposited over a time ~ 1 s after core. (Likewise, 
there will be a corresponding, though far smaller increase 
in negative mu and/or tau lepton number stemming from 

The I'e's deposited by i^e could represent a signifl- 
cant increase in electron lepton number. In turn, this will 
tend to decrease the neutron excess. This is because the 
additional t/g's will tend to shift the equilibrium relation, 
Ue + n ^ p + , to the right, producing higher electron 
fraction Y^ and more protons. Likewise, the neutron- 
to-proton ratio, n/p = Y^^ — 1, in the material near 
the neutrino sphere will be transmitted by the and 
fluxes emergent from the neutrino sphere to the material 
in the region between the neutron star and the shock [soj . 

In the early shock re-heating regime, this increase in 
electron fraction Ye in the material ejected by neutrino 
heating could be beneficial for nucleosynthesis. In the 
calculations of nucleosynthesis in early, shock re-heating 
epoch neutrino-heated ejecta performed by Woosley et al. 
using the Mayle and Wilson supernova simulation results 
[6d |, it was found that there was an overproduction of 
neutron number TV = 50 nuclei. Subsequently it was 



pointed out in Ref. [6l| that a modest increase in Yg could 
cure this problem. The Vg lepton number 

transfer process at least sends Yg in the right direction 
at the right epoch to help. 

Effects of active-sterile-active neutrino flavor transfor- 
mation at later times, in the post-shock revival hot bub- 
ble, may be very interesting, but are beyond the scope of 
the current work. 

We note that re-conversion of sterile neutrinos has 
been considered pre viously in models for r-proccss nu- 
cleosynthesis [62l. l63l| ■ These calculations, however, con- 
centrated on the late-time regime above the core and con- 
sidered a much different sterile neutrino mass and vac- 
uum mixing range from the one considered here. Addi- 
tionally, active-active neutrino flavor transformation in 



the supernova environment is a very difficult problem 

[H m, m, [el [H, [13, H [H [zO, [zi [zi [zl, Izl, l^. 

A complete assessment of nucleosynthesis effects would 
necessitate treating all active-active and active-sterile 
neutrino flavor conversion processes. 



VI. CONCLUSIONS 

Generally it has been assumed that the emission of 
sterile neutrinos from the supernova core will tend to de- 
crease the prospects for obtaining a successful core col- 
lapse supernova explosion. This may be true if a large 
enough amount of energy is lost from the core. This is 
because, after all, most of the gravitational binding en- 
ergy released in the collapse of the core and subsequent 
quasi-static contraction of the hot proto-neutron star is 
"stored" in trapped seas of active neutrinos of all species. 
Moreover, it is this neutrino energy which, ultimately, 
will be invoked one way or another to revive the nuclear 
photo-dissociation-degraded bounce shock. 

However, in this paper we point out that the notion 
that sterile neutrino emission is bad for shock revival is 
predicated on the assumption that there will be no re- 
conversion of these sterile neutrinos to active neutrino 
species. Indeed, our calculations suggest that such a re- 
conversion process could take place under some circum- 
stances and that this re-conversion could effect an en- 
hancement in energy and electron lepton number trans- 
port from deep in the core to the regions just below the 
neutrino sphere. This could increase the prospects for a 
viable explosion through: (1) pre-heating of the material 
ahead of the shock causing a reduction in the nuclear 
photo-disintegration burden on the shock; and (2) en- 
hancement of the Ve and v^, heating rate of the material 
under the bounce shock. 

We have found that the sterile neutrino mass and mix- 
ing parameters for which these enhancement processes 
can take place conform to our earlier estimates [32] of 
these: sterile neutrino rest mass range 1 keV < ^ 
5 keV; and v^. ^ effective 2x2 vacuum mixing angle 
in a range satisfying 10~^° < sin^ 26* < 10~^. Most sig- 
nificantly, we find that the neutrino mass and mixing pa- 
rameter ranges which give supernova explosion enhance- 
ment include those ranges of parameters which give a 
possibility for viable sterile neutrino dark matter. What 
was missing in our earlier work i33] was an assessment 
of the effects of the shock itself on the neutrino forward 
scattering potential which governs active-sterile neutrino 
flavor transformation. In this paper we have done this 
assessment. 

However, there are many uncertainties and our one- 
zone calculations can be regarded only as rough outlines 
for how active-sterile-active neutrino flavor conversion 
processes affect supernova core and shock physics. How 
can our calculations be improved on? 

First, in the context of a realistic proto-neutron star 
model, a self consistent hydrodynamic treatment of shock 
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propagation coupled with active-sterile and sterile-active 
neutrino flavor transformation processes is in order. This 
could resolve tricky issues associated with the effec- 
tiveness of pre-heating in relieving the nuclear photo- 
dissociation burden on the shock. 

Second, it would be useful to employ a detailed treat- 
ment of neutrino transport, coupled with a realistic 
model for the structure and equation of state of the region 
of the proto-neutron star near the neutrino sphere, to as- 
sess the way in which energy deposited via u^, ^ i^g ^ i^e 
is divided up among the various active neutrino species. 
Also, we need to know how this deposited energy affects 
Ye and the emergent luminosities of the active neutrino 
species at and above the neutrino sphere. 

There is yet a third source of uncertainty, one which 
may be an issue for all core collapse supernova models. 
We have pointed out in this paper that the initial core 
bounce shock strength is an important quantity for char- 
acterizing how the shock modifies the "fossil" neutrino 
forward scattering potential profile which is left at the 
end of the core in-fall epoch. The initial shock strength 



depends on many factors in both the pre-coUapse hydro- 
static evolution epochs of the progenitor star as well as on 
in-fall physics issues like nuclear weak interaction rates 
and the sub-nuclear density equation of state. 

Ultimately, of course, the core collapse supernova prob- 
lem is a grossly nonlinear one. We will have to grapple 
with this nonlinearity, as well as a host of fundamental 
nuclear physics and multi-dimensional hydrodynamic is- 
sues, if we ever hope to realize the awesome power of 
this "laboratory" for revealing/constraining new physics 
beyond the Standard Model. 
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